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We analyze the sensitivity of gravitational-wave antenna with stable double optical spring cre- 
ated by two independent pumps. We investigate regime of three close eigen frequencies (roots of 
characteristic equation) which appears to provide more wide frequency band in which sensitivity of 
antenna can beat Standard Quantum Limit (SQL) than previously considered regime with two close 
eigen frequencies. We take into account optical losses and show that they do not degrade sensitivity 
significantly. We also demonstrate possible application of considered regime to Einstein Telescope. 



I. INTRODUCTION 

Currently the search for gravitational radiation from 
astrophysical sources is carried out with the first- 
generation Earth-based laser interferometers (LIGO in 
USA [1-3], VIRGO in Italy [4, 5], GEO-600 in Germany 
[6, 7], TAMA-300 in Japan [8, 9] and ACIGA in Australia 
[10, 11]). The development of the second-generation GW 
detectors (Advanced LIGO [12, 13], Advanced Virgo [14], 
GEO-HF [15] and LCGT [16]) is well underway. 

The sensitivity of the first-generation detectors is lim- 
ited by noises sources of various nature: seismic and sus- 
pension thermal noise at low frequencies (below ~ 50 
Hz), thermal noise in suspensions, bulks and coatings of 
the mirrors (~ 50 — 200 Hz), photon shot noise (above 
~ 200 Hz). It is expected that the sensitivity of the 
second-generation detectors will be ultimately limited by 
the noise of quantum nature arising due to Heisenberg's 
uncertainty principle [17-20] over most of the frequency 
range of interest. The optimum between measurement 
noise (photon shot noise) and back-action noise (radi- 
ation pressure noise) is called the Standard Quantum 
Limit (SQL) This level is expected to be reached in 
the forthcoming second generation of large-scale laser- 
interferometric gravitational-wave detectors. Third gen- 
eration detectors, such as the Einstein Telescope [21] aim 
to significantly surpass the SQL over a wide frequency 
range [22]. 

The most promising methods to overcome the SQL rely 
on the implementation of optical (ponderomotive) rigid- 
ity [20, 23-25] which effectively turns the test masses of 
a gravitational-wave detector into harmonic oscillators 
producing gain in sensitivity [26-31]. A single optical 
spring always causes instability. The instability (nega- 
tive damping for mechanical degree of freedom) can be 
compensated by incorporating a linear feedback control 
loop and in the ideal case (no additional noise is intro- 
duced by the feedback) it would not modify the noise 
spectrum of a GW detector [29]. In practice, however, 
the need for control gain at frequencies inside the detec- 
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FIG. 1: Scheme of an Advanced LIGO interferometer pumped 
by two lasers. The main laser is detuned to give a positive 
optical rigidity and negative damping (it is tuned on the right 
slope of resonance curve) while the auxiliary laser is detuned 
to give negative optical rigidity and positive damping. 



tion band can cause undesirable complexity in the control 
system or can introduce additional classical noise. 

An alternative way to suppress the instability was pro- 
posed [32] and experimentally demonstrated [33], by in- 
jecting a second carrier field from the bright port (see 
Fig. 1) in order to create a relatively small additional neg- 
ative rigidity component, thus leaving the total rigidity 
(of both lasers together) positive, but at the same time to 
introduce a relatively large additional positive damping 
component to make total damping positive. The main 
purpose of the second carrier is to create a second op- 
tical spring that forms a stable optical spring together 
with the first one — even though each individual optical 
spring, acting alone, would be unstable. Both carriers are 
assumed to have different polarizations (or carriers' fre- 
quencies to differ from each other by large enough value) , 
so that there is no direct coupling between the two pump 
field (although they both directly couple to the mirrors) . 

A simple criterion for the stability of an optical spring 
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and its application to the double resonance regime was 
presented in [34]. It allows to use attractive regimes of 
double resonance [30, 31, 34] or negative inertia [35] in 
stable variants with no need to use feedback loops. 

In this paper we further analyze the sensitivity of 
Advanced LIGO (aLIGO) interferometer with a double 
pump using approach [32] . To obtain more wide sensitiv- 
ity curve we use the regime of three close roots of char- 
acteristic equation instead of double resonance regime. 
However this regime demands smaller interferometer re- 
laxation rates, close to limit caused by optical losses, 
so there is need to consider optical losses too. Analy- 
sis shows that presence of losses doesn't affect the sensi- 
tivity curves much. Finally we apply this regime to the 
parameters of Einstein Telescope [22] and show that use 
of stable double optical spring makes possible to make an 
improvement to sensitivity in small frequency range. 
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II. OUTPUT ANALYSIS 

We consider the balanced interferometer similar to 
aLIGO configuration (see Fig. 1) pumped by two lasers. 
We assume that vacuum fluctuations input through dark 
port (no squeezing), Fabry-Perot cavities in arms are 
identical and have no optical losses. The masses of input 
and end mirrors in arms are equal to m. The mean fre- 
quencies of each pump w\ t 2 are equal to one of the eigen 
frequencies of each cavity. Here and below subscripts 12 
refer to main and auxiliary pumps. We also introduce 
the notations: 
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Here L is the distance between the mirrors in arms (4 km 
for aLIGO), £ is the distance (several meters) between SR 
mirror and the input mirrors in Fabry-Perot cavities, due 
to strong inequality I -C L we assume that 4>i,2 do not 
depend on VL. We assume that arm cavities are tuned in 
resonance, 70 is the relaxation rate of single Fabry-Perot 
cavity in arm, T is amplitude transmittance of input mir- 
rors in arms, p is amplitude reflectivity of SR mirror, Ai t 2 
are the detunings introduced by displacement shift of SR 
mirror, 1^2 is the relaxation rates of the differential mode 
of the interferometer for each pump. It is worth under- 
lying that detunings A1.2 and relaxation rates Ti^ may 
be different for each pump. 

We start from equations for output quadratures am- 
plitudes 6 1 1 2 2 ' ) expressed in terms of input quadratures 
amplitudes cl^^ m frequency domain (superscripts < - 1,2 ) 



Here \1/ is susceptibility of mechanical degree of freedom 
accounting optical rigidities K\ 2 introduced by pumps 1 
and 2, Ii 2 are mean powers circulating in arms (pumped 
by main and auxiliary lasers), p is reduced mass, n 
is spectral frequency, h is dimcnsionless gravitational 
metric perturbation normalized by SQL perturbation 
/isql- The term J- in formulas (2) describes fluctuational 
(back action) and signal forces, while the first terms in 
right parts describe measurement errors. The coefficients 
Ui, Tj are the following: 
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The equations (2,3) slightly differ from equations used in 
[32] due to different definition of quadrature amplitudes 
(see (Al, A2)). 

In experiment homodyne detector of output field mea- 
sures arbitrary combination of quadratures in each chan- 
nel defined by homodyne angles £i i2 : 
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For input field in vacuum state quadratures a\ do 
not correlate to each other and their single-sided spectral 
densities [36] are equal to 
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In order to find condition under which sensitivity is 
better than SQL one may apply the following semiquali- 
tative consideration. 

From equations (2) we see that if susceptibility is large 
enough (|/if2 2, L| 1) the back action force prevails over 
measurement errors and mainly defines sensitivity and 
spectral density ShBA of noise recalculated to dimension- 
less metric h is (see Eq. 2e): 
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It is obvious to assume that close to resonance (when 
l^il 2 ^! 3> 1) we have approximation: \Ki + K 2 \//j,n 2 ~ 
1. Hence, from (6) one may conclude that this term will 
be less than unity (i.e. sensitivity is better than SQL, 
S'hBA "C ^sql) m case °f small optical relaxation rates, 
i.e. 
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If rigidity |_ffi| introduced by main pump is larger than 
\K 2 \ we should keep in mind that first inequality (Fi <C 
Ai) in (7) is more important than second one (r 2 <C A 2 ). 



TABLE I: Planned parameters of Advanced LIGO 



Parameter 


Value 


Arm length, L 


4 km 


Mass of each mirror, m 


40 kg 


ITM amplitude transmittance, T 


V0.005 (V0.015) 


Relaxation rate of single FP cavity, 70 (la) 


94 s _1 


SRM amplitude reflectivity p 


VoM (Vols) 


Optical wavelength, A 


1064 nm 


Optical losses A 2 in each arm per roundtrip 


10" 5 



For frequencies far from resonance susceptibility de- 
creases and when |/if2 2, I'| ~ A1/T1 back action and 
measurement noises will be approximately equal to each 
other. Hence, one may get better than SQL sensitivity 
in bandwidth where susceptibility is high enough. 

In more general case formula for sensitivity is more 
complicated than (6). If we measure, for example, only 



quadrature the sensitivity is equal to 



S. 



(C) 



h 2 

SQL 



cos£ 



nD, 2 ^f(Ui cosC + U 2 sin C) 



sin( 



/ifi 2 *(C/iCosC + C/ 2 sinC) 



•T 2 



(8) 



|T 3 | 2 T 4 | 2 
2 2' 



2.00 

1.50 
1 .00 
0.70 
0.50 

0.30 

0.20 
0.15 




FIG. 2: Trace 1: sensitivity £ c -° (/) (8); 2: sensitivity £ba(/) 
(6) defined only by back action; 3: corresponding dimension- 
less susceptibility tp = fi(2nf) 2 ^/ for parameters (9). 



Note that for single-pumped interferometer regime with 
small relaxation rate is not the attractive one because it 
may provide SQL overcoming only in small bandwidth 
(decreasing with decrease of Ti). That is why in aLIGO 
relatively small transmittance yl — p 2 is planned. From 
formulas (lc, Id) and parameters listed in Table I one 
may easy estimate that ratio Ti/Ai is restricted by 
Ti/Ai > 0.025. Below we do not take into account this 
restriction assuming that parameters Fi. 2 , Ai. 2 may be 
chosen arbitrary and transmittances T and \]\ — p 2 may 
be changed. 



A. Close resonances regime 

In double-pumped interferometer one may manipulate 
by value of susceptibility "J in wide range due to sophisti- 
cated frequency dependence of optical string. It is impor- 
tant that this could be done avoiding instability (it means 
that imaginary parts of roots of characteristic equation 
= should have negative signs) [34]. One of possi- 
bilities is the usage of so called double resonance regime 
[27, 30, 31] which allows to obtain two close or coinciding 
resonance frequencies. 

In order to obtain more wide bandwidth of high sus- 
ceptibility one may use the case of three resonance fre- 
quencies being relatively close to each other. As example 
we present the plots of sensitivity and susceptibility on 
Fig. 2 for such particular case. For these plots we used 
the following parameters 



7i~731kW, ri~3.06s" 
7 2 ~1.06kW, r 2 ~15.3s~ 



Ai ~ 961s- 1 , 
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For these parameters the complex roots of characteristic 
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equation are equal to 

fti ~ ±2vr x 96.6 - i 0.49 s" 1 , (10a) 
Vt 2 ^ ±2tt x 102.4- i 0.63 s" 1 , (10b) 
n 3 ~ ±2vr x 116.3- il. 81s" 1 . (10c) 

The real part of roots are eigen frequencies and imagi- 
nary parts describe relaxation (negative sign of imaginary 
parts corresponds to relaxation, positive one — to insta- 
bility). All roots (10) are stable. We used procedure of 
root calculation [34] allowing to have stable roots (with 
negative imaginary parts) avoiding instability. 

Strictly speaking we should use both outputs measur- 
ing the quadratures (4) and taking their weighted sum 
(the corresponding procedure is given in Appendix to 
[32]). However, for our case (9) the circulating powers 
Ji and I2 differ to each other so dramatically that useful 
information in port 2 is negligible small and it is quite 
enough to measure only quadrature b\. In particular, on 
Fig. 2 we present sensitivity (8) for the case C, = 0. 

The presented example demonstrates the main advan- 
tage of close resonances regime: high sensitivity (about 
10 times better than SQL) in bandwidth as wide as about 
one half of mean frequency. 

Note that sensitivity curve presented on Fig. 2 may be 
easy tuned to another frequency range by control powers 
/i,2, detunings Ai^ and relaxation rates T\p. For ex- 
ample, for frequency range around 30 Hz one may easy 
recalculate: 




FIG. 3: Traces 1,2: sensitivities f^ =0 (/) without and with ac- 
count of optical losses (Fa = and Ta = T\ correspondingly); 
3: sensitivity £ba (6) defined only by back action; 4: dimen- 
sionless susceptibility tp = /i(27r/) 2 ^I' for parameters (11). 



where Ti describes relaxation through SRM (lc) (may 
be called as load relaxation) and Ta — relaxation due to 
optical losses (loss relaxation). 

There is a useful rule [31]: a lossy optical position me- 
ter is equivalent to the similar lossless one with gray filter 
(with power transmission Ta/Tx) attached to its signal 
port. It means that for lossy case, for example, the for- 
mula (2) for output quadrature has to be rewritten 
in form: 
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This case is more attractive due to modest requirements 
for powers circulating in arms (18 kW instead of 730 kW). 
However, estimates above shows tough requirements for 
bandwidth Fx ~ 0.89 s" 1 , which is close to limit due to 
inevitable optical losses. This fact forces us to consider 
optical losses. 



1. Account of optical losses 

We will take into account only the optical losses in 
arms' mirrors which may be characterized by round trip 
loss coefficient A 2 . We do not take into account losses in 
beam splitter and SRM because its' influences are negli- 
gible small. 

It is known [31] that relaxation rates Tip of differen- 
tial mode of interferometer with optical losses (shown on 
Fig. 1) may be separated into two parts 
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where quadrature describes additional vacuum fluc- 
tuations (with the same spectral density (5)) appearing 
due to optical losses. Generalization for other quadra- 
tures is obvious. Additionally we should replace in all 
formulas for 'J, Ui, Ti relaxation rates Ti^ by Ti^- 

For planned in aLIGO losses (see Table I) one may 
estimate Ta — 0.4 s" 1 . In according with estimate (11) 
we have to have r\ ~ 0.89 s" 1 , i.e. the load and loss 
relaxations are approximately equal: Ti ~ Ta- 

The corresponding sensitivities plots are presented on 
Fig. 3. We see that so dramatically large optical losses 
practically do not bring any significant degradation to 
the sensitivity (!). 

So we see that it is possible to circumvent SQL by 
about 10 times in wide bandwidth (about the mean fre- 
quency) using two optical stable springs. 

Note that for simplicity we considered above not op- 
timal case when only quadrature in one channel is mea- 
sured. However, this sensitivity only slightly differs from 
maximal sensitivity when quadratures from both chan- 
nels with optimal homodyne angles are measured and 
then are properly combined — see details in Appendix B. 
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if the following inequality is fulfilled 
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FIG. 4: Plots of sensitivity £ of ET normalized to SQL ver- 
sus frequency. Traces 1,2: different regimes of optical rigidity. 
Trace 3: design sensitivity of ET without squeezing, no ther- 
mal noise present. 



2. Application to Einstein Telescope 

As the most promising device for gravitational wave 
detection in future is going to be Einstein Tele- 
scope (ET), it seems useful to try out double optical 
spring for its configuration. 

We take the parameters for low-frequency setup of ET 
(see [22]) and recalculate regime from fig. 3 using these 
parameters and aiming to obtain mean frequency approx- 
imately equal to 7 Hz. The result of this recalculation 
is trace 1 in fig. 4. For each of traces 1 and 2 at fig. 4 
we define optimal set of quadratures and filter functions 
to maximize sensitivity according to algorithm given in 
appendix B. 

The regime of three close roots however occurs to pro- 
vide very narrow sensitivity curve as compared with ET 
design configuration (trace 3 at fig. 4). Using procedure 
described in [34] it is quite easy to find some interme- 
diate regimes which demonstrate wider sensitivity curve 
than trace 1 does, but with peak sensitivity better then 
one designed for ET. We present one such regime plot- 
ted by means of trace 2 at fig. 4. These regimes mainly 
differ by detunings and relaxation rates of main and aux- 
iliary pumps. Required values of power in main pump are 
approximately equal for traces 1 and 2. Also these val- 
ues are smaller than power in design configuration (18 
kW). However power of auxiliary pump changes signif- 
icantly (still remaining small as compared to power of 
main pump) while sensitivity curve becomes wider. For 
details see Appendix C. 



B. Speculations on quantum behavior of 
optomechanical degrees of freedom 

Let remind that oscillator in thermal bath at temper- 
ature T will demonstrate quantum behavior [17, 20, 37] 



< 1, 



(14) 



where h is Boltzmann constant, Q and Q arc frequency 
and quality factor of oscillator. 

The plots of susceptibility on Figs. 2, 3 show three 
peaks corresponding to three optomechanical degrees of 
freedom. It is important that they are stable ones. Gen- 
erally speaking it is not clear how to describe correctly 
these degrees of freedom similar to normal coordinates 
in case of coupled oscillators. However, for estimation 
one can consider the mechanical degree of freedom and 
calculate mean kinetic and potential energy Ei stored in 
each peak: 
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where Qi and AQi are mean frequency and bandwidth of 
each peak (real and imaginary parts of roots (10,11)). Let 
assume that only fluctuational light pressure forces orig- 
inate fluctuations of displacement z (no thermal noise). 
The numerical calculations gives the following values: 

~ 1.28ftf2 1; E 2 ~1.76hQ 2 , J5 3 ~l.lfin 3 (16) 

We see that mean energy in each degree of freedom is 
close to one quanta, hence, they should demonstrate 
quantum behavior. It becomes obvious if we apply crite- 
ria (14) for each peak: 
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(17) 
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It provides incredible possibility to observe quantum 
behavior of really macroscopic object with mass y, ~ 
10 kg (!) 



III. CONCLUSION 

We have shown that usage of stable double optical 
spring regime in laser gravitational detectors may pro- 
vide sensitivity gain if relaxation rates of optical modes 
are much smaller than detunings. However, the decrease 
of relaxation rates is restricted by optical losses in mir- 
ror. We have shown that even for case when relaxation 
through signal recycling mirror (load relaxation) and re- 
laxation via optical losses (loss relaxation) are equal to 
each other the degradation of sensitivity is relatively 
small. 

The presented example of close resonances regime via 
double pumped optical spring promises the possibility to 
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circumvent Standard Quantum Limit by about ten times 
in the frequency range about half of mean frequency. 
Experimenter may further vary parameters set (pump 
powers, detunings) to control susceptibility and, hence, 
sensitivity curve. Of course, this gain of sensitivity will 
take place if the level of thermal and technical noises is 
low enough. The thermal noise in mirror's coating makes 
the main contribution in thermal noise budget. However, 
one may hope for progress in manufacture of the interfer- 
ometric coating. Note that the planned level of thermal 
noise in laser gravitational detectors of third generation 
(Einstein Telescope ET-D Low Frequency Interferome- 
ter) is planned to be about ten times less than Standard 
Quantum Limit [38] in frequency range about 10 Hz. 

If thermal noise is small enough the stable double opti- 
cal spring provides unique possibility to observe quantum 
behavior of macroscopic object with effective mass about 
10 kg, because fluctuational energy (created by back ac- 
tion force) stored in each peak of sensitivity is about one 
quanta Mli. 
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Appendix A: Initial formulas 

In this Appendix we derive formulas (2, 3). The elec- 
tric fields E\ >2 in propagating wave of each pump and 
corresponding the mean intensities J\^ 2 of light beam can 
be written as follows [361: 



2tt fiwi : 2 
Sc 



E li2 ~ A/ -';"^ e-^x 



x Ai 2 + 



Ji,2 = t^iMia\ 2 , [<H,2(n), a+ 2 (Q')] = 2n5(n - Q'), 

where S is the cross section of the light beam, c is the 
velocity of light, a and a + are annihilation and creation 
operators. 

One can obtain the formula in frequency domain for 
the output fields 61.2 in dark port as a function of the 
input fluctuational field ai <2 and mirror positions [32]: 
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where xe, Ve, Xn, Vn — are the displacements of FP 
mirrors (see notations on Fig.l), I± t2 are mean optical 
power circulating in each arm from pumps 1 or 2. 

As we assume that input fields are in vacuum state (no 
squeezing), it is convenient to denote 
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As usual we introduce input a\ 2 and output b\ 2 
quadrature amplitudes 

a$ = [5i,a(fi) + af t3 (-fl)] /V2, (A3a) 

afl = [ai l2 (fi) - a+ 2 {-Sl)] /iV2, (A3b) 

b^ 2 = \b 1 , 2 (Q) + b+ 2 (-n)]/V2, (A3c) 

b<?l = 161,2(0) - b+2(~ n )] A^- ( A3d ) 

Now we have to account evolution of differential coor- 
dinate z through susceptibility ^ in frequency domain: 

z=^(Q)(h + f 2 +F s ), F s =^ 2 Lh. (A4) 

Here F s is the equivalent signal force. Back action forces 
/i,2 produced by fluctuations of light pressure are equal 
to: 
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Substituting (Al, A2, A4, A5) into definition (A3) of 
quadrature amplitudes one may obtain the formulas (2, 
3). 



Appendix B: Accurate sensitivity calculation 

The sensitivity of GW-antenna depends on homodyne 
angles and the filter functions that are used during pro- 
cession of quadratures and homodyne currents. The op- 
timal set of filter functions can be found as the eigen 
vector of the proper matrix [32, 39]. Below we provide 
the algorithm to obtain this set in general case of multi- 
ple pumps and calculate optimal sensitivity for the case 
of two pumps. 

Consider that measurement of output quadratures 
gives several homodyne currents ji (defined by formu- 
lae (4) in case of two pumps) , each of which consists of 
signal part and noise part j> N '. The total output 
is built as the weighted sum of these currents with filter 
functions Yr. 
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FIG. 5: Plots of sensitivity versus frequency. 1: optimized 
sensitivity £,(/), losses are absent; 2: sensitivity £ °(/) in 
case of measurement of only the amplitude quadrature of the 
strong pump, losses are absent. 3 and 4: the same with losses 

r A = r 1 . 

Spectral density of signal S s and noise S n parts of output 
signal can be presented as the quadratic forms over the 
vector of filter functions Y = {Yi; Y 2 ; . . . }. 

S n = (Y;NY) = J2ytN ij Y j ; 

ij 

S S = (Y;SY)^J2 Y * S ^ 

The elements of matrices S and N of the corresponding 
quadratic forms are defined by expressions 



WCViPif; N ij{ n,0=r>jf\ (Bi) 

where the line over means the operation of calculation 
of the cross spectral density. These elements depend on 
frequency il and homodyne angles used in quadratures 
measurements . 

Inverse sensitivity of antenna is given as the ratio 

"•SQL °n 

It depends on the same quantities and also on the set 
of filter functions. The optimal set of the latter should 
maximize thus minimizing the sensitivity. 

The eigen vector of the matrix P = N^ 1 ■ S that cor- 
responds to it's largest eigen value provides desired set 



of filter functions [32, 39]. Corresponding sensitivity is 
given by inverse to that eigen value. 

The optimal sensitivity calculated following this 
method depends on frequency and homodyne angles. Op- 
timization (minimization) over the latter gives optimal 
sensitivity for the case of frequency dependent ho- 

modyne angles. Due to experimental difficulties arising 
of realisation of such regime it is worth to make estima- 
tion for optimal set of homodyne angles. 

As a criterion we use the condition of minimization of 
integral sensitivity S 

n 2 

5(c) = J m()dn, 

where frequencies Oi ; 2 are boundaries of the range in 
which the sensitivity £^ =0 (r2) is under the SQL level. 
This optimization gives the best sensitivity £j acheivable 
with frequency independent homodyne angles. 

In the case of double pump all the matrices and vectors 
are two-dimensional so all the eigen vectors and eigen val- 
ues can be calculated analytically. Further optimization 
over the homodyne angles is done numerically. Results 
are presented in the fig. 5. One can see that that sensi- 
tivity does not differ dramatically as compared with the 
approximation presented in Sec. II (taking into account 
only output from one pump). 



Appendix C: Tuning details for fig. 4 

Parameters needed to achieve regimes with sensitivi- 
ties represented by first two traces at fig. 4 are listed in 
Table II below. 

TABLE II: Parameters for plots at fig. 4 





Trace 1 


Trace 2 


h, kW 


4.5 


5.4 


h, kW 


0.06 


0.12 


At, s" 1 


65.7 


70.5 


A 2 , s" 1 


-42.8 


-43.8 


Pi, s" 1 


0.21 


1.05 


r a , s- 1 


1.04 


6.70 
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